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Abstract 

We investigate the vacuum energy density induced by quantum fluctuations of a bulk scalar 
field with general curvature coupling parameter on two codimension one parallel branes in a 
(D + l)-dimensional background spacetime AdSm+i x £ with a warped internal space S. It 
is assumed that on the branes the field obeys Robin boundary conditions. Using the general- 
ized zeta function technique in combination with contour integral representations, the surface 
energies on the branes are presented in the form of the sums of single brane and second brane 
induced parts. For the geometry of a single brane both regions, on the left (L-region) and on the 
right (R-region), of the brane are considered. The surface densities for separate L- and R-regions 
contain pole and finite contributions. For an infinitely thin brane taking these regions together, 
in odd spatial dimensions the pole parts cancel and the total surface energy is finite. The parts 
in the surface densities generated by the presence of the second brane are finite for all nonzero 
values of the interbrane separation. The contribution of the Kaluza-Klein modes along E is in- 
vestigated in various limiting cases. It is shown that for large distances between the branes the 
induced surface densities give rise to an exponentially suppressed cosmological constant on the 
brane. In the higher dimensional generalization of the Randall-Sundrum braneworld model, for 
the interbrane distances solving the hierarchy problem, the cosmological constant generated on 
the visible brane is of the right order of magnitude with the value suggested by the cosmological 
observations. 
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1 Introduction 

Motivated by the problems of the radion stabilization and the generation of cosmological constant, 
the role of quantum effects in braneworlds has attracted great deal of attention pQ-[47]. A class of 
higher dimensional models with the topology AdS/^+i x S, where £ is a one-parameter compact 
manifold, and with two branes of codimension one located at the orbifold fixed points, is considered 
in Refs. [23, 24J. In both cases of the warped and unwarped internal manifold, the quantum effective 
potential induced by bulk scalar fields is evaluated and it has been shown that this potential 
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can stabilize the hierarchy between the Planck and electroweak scales without fine tuning. In 
addition to the effective potential, the investigation of local physical characteristics in these models 
is of considerable interest. Local quantities contain more information on the vacuum fluctuations 
than the global ones and play an important role in modelling a self-consistent dynamics involving 
the gravitational field. In the previous papers \40\ [4"T] we have studied the Wightman function, 
the vacuum expectation values of the field square and the energy-momentum tensor for a scalar 
field with an arbitrary curvature coupling parameter obeying Robin boundary conditions on two 
codimension one parallel branes embedded in the background spacetime AdS^+i x S with a warped 
internal space S. For an arbitrary internal space S, the application of the generalized Abel-Plana 
formula [48j allowed us to extract form the vacuum expectation values the part due to the bulk 
without branes and to present the brane induced parts in terms of exponentially convergent integrals 
for the points away from the branes. 

The braneworld corresponds to a manifold with boundaries and the physical quantities, in 
general, receive both volume and surface contributions. In particular, the contributions located on 
the visible brane are of special interest as they are direct observables in the theory. In Ref. [32] 
the vacuum expectation value of the surface energy-momentum tensor is evaluated for a massive 
scalar field subject to Robin boundary conditions on two parallel branes in (D + l)-dimensional 
AdS bulk. It has been shown that for large distances between the branes the induced surface 
densities give rise to an exponentially suppressed cosmological constant on the brane. In the 
Randall-Sundrum braneworld model, for the interbrane distances solving the hierarchy problem 
between the gravitational and electroweak mass scales, the cosmological constant generated on 
the visible brane is of the right order of magnitude with the value suggested by the cosmological 
observations. In the present talk based on |47j . we describe similar issues within the framework of 
higher dimensional braneworld models with warped internal spaces. 

2 Surface energy-momentum tensor and the zeta function 

Consider a scalar field <p(x) with curvature coupling parameter £ satisfying the equation of motion 

{V M V M + m 2 + (R) <p{x) = 0, (1) 

where Vm is the covariant derivative operator and R is the scalar curvature for a (D+l)-dimensional 
background spacetime. We will assume that the bulk has the topology AdS^+i x S and is described 
by the line element 

ds 2 = g MN dx M dx N = e~ 2kDy U^daPdaf - j^dX^X^ - dy 2 , (2) 

where rj^ is the metric tensor for Di-dimensional Minkowski spacetime, hp is the inverse AdS 
radius, and the coordinates X 1 , i = 1, . . . , D2, cover the internal manifold E, D = D\ + Z?2- The 
Ricci scalar corresponding to line element ([2]) is given by formula R = —D(D + l)k 2 D - e 2kD yR (j) , 
with Rm being the scalar curvature for the metric tensor 7^. In the discussion below, in addition 
to the radial coordinate y we will also use the coordinate z = e k ° y /&£>, in terms of which the line 
element is written in the form conformally related to the metric in the direct product spacetime 
R( D t< 1 ) x E by the conformal factor (kr>z)~ 2 . 

We are interested in one-loop vacuum effects induced by quantum fluctuations of the bulk field 
(f(x) on two parallel branes of codimension one, located at y = a and y = b, a < b. We assume 
that on the branes the field obeys Robin boundary conditions 

[Ay + Bydy^j Cf(x) =0, ?/ = O, 6, (3) 
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with constant coefficients A y , B y . In a higher dimensional generalization of the Randall-Sundrum 
braneworld model the coordinate y is compactified on an orbifold S 1 /Z2, of length I, —I ^ y ^ I, 
and the orbifold fixed points y = and y = I are the locations of two branes. The corresponding 
line-element has the form ([2]) with the warp factor e~ 2ko ^. In these models the region between 
the branes is employed only. For an untwisted bulk scalar with brane mass terms c a and q,, the 
corresponding ratio of the coefficients in the boundary condition ([3]) is determined by the expression 
(see, e.g., Refs. [50l El [35] for the case of the bulk AdS/j+i and Refs. [23l [40] for the geometry 
under consideration) 

Aj n^cj + 4D(k D 
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(°) = l, n (b) = -1. (4) 



In the supersymmetric version of the model [50j one has q> = — c a and the boundary conditions 
are the same for both branes. For a twisted scalar, Dirichlet boundary conditions are obtained on 
both branes. 

For the geometry of two parallel branes in AdS^+i x E with boundary conditions ([3j), the 
Wightman function and the vacuum expectation values (VEVs) of the field square and the bulk 
energy-momentum tensor are investigated in Refs. [40, 41J. On manifolds with boundaries the 
energy-momentum tensor in addition to the bulk part contains a contribution located on the 
boundary. For an arbitrary smooth boundary dM with the inward-pointing unit normal vector 
n L , the surface part of the energy- momentum tensor for a scalar field [51] is given by the formula 

^mn = $( x 'i 9M)tmn j where the "one-sided" delta-function 5(x;dM) locates this tensor on dM 
and 

tmn = (ip 2 K M N ~ (2C - l/2)h M N¥n L X7 L (p. (5) 

In Eq. ©, hj\jN = 9mn + ^M^jv is the induced metric on the boundary and Kmn = hj^h^V Lnp 
is the corresponding extrinsic curvature tensor. 

By expanding the field operator over a complete set of eigenfunctions {ip a (x), ¥* a {x)} obeying the 
boundary conditions and using the standard commutation rules, for the VEV of the operator tmn 
one finds (0\tmn\0) = J2 a T MN{ l Pa(x), </?*(x)}, where the bilinear form tmnW,^} is determined 
by the classical energy- momentum tensor ([5]). Below in this section, we will consider the region 
between the branes. The corresponding quantities for the regions y ^ a and y ^ b are obtained as 
limiting cases. As we have mentioned before, in the orbifolded version of the model, the only bulk 
is the one between the branes. In this region the inward-pointing normal to the brane at y = j, 
j = a, b, and the corresponding extrinsic curvature tensor are given by the relations n^ M = n&5% 

and K^ N = —n^kahMN, where is defined in formula ([3]). By using these relations and the 
boundary conditions, the VEV of the surface energy-momentum tensor on the brane at y = j is 
presented in the form 

(0|r[^|0) = -h M Nn^k D C^\ip 2 \Q,) y=j , C s = (- (2C - l/2)Aj/(kr>Bj). (6) 

From the point of view of physics on the brane, Eq. ([6]) corresponds to the gravitational source 
of the cosmological constant type with the surface energy density = (0|tq^°|0) (surface energy 
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the equation of state ej v = —pf- It is remarkable that this relation takes place for both subspaces 



per unit physical volume on the brane at y = j or brane tension), stress = — (0|r} J |0), and 
-W _ » 



on the brane. It can be seen that this result is valid also for the general metric instead of rj^ 
in line element ([2]). For an untwisted bulk scalar in the higher dimensional generalization of the 
Randall-Sundrum braneworld based on the bulk AdS^+i x S the coefficient in Eq. ([6]) is given by 
the formula 

C 3 = C(4DC - D + 1) + (C - l/4)nW Cj /k D , (7) 



3 



In particular, the corresponding surface energy vanishes for minimally and conformally coupled 
scalar fields with zero brane mass terms. 

In order to evaluate the expectation value of the surface energy-momentum tensor we need the 
corresponding eigenfunctions. These functions can be taken in the decomposed form 

Va {x M ) = f ?} y)e ~T\ k * MX), & = ("A»,k), W/9,„ = Jk* + ml + \% k = |k|, (8) 

where the mass spectrum m n is determined by the boundary conditions. The modes ipp{X) are 
eigenfunctions for the internal subspace: 

[A (7) + 0R (7) ] MX) = -A^(X), (9) 

with the eigenvalues A| and the orthonormalization condition J dP 2 X y i /jij;p(X)ip%,(X) = 8ppi. In 
Eq. (jUJ), ^(7) is the Laplace-Beltrami operator for the metric 7^. Prom the field equation JT]), one 
obtains the equation for the function f n {y) with the solution 

Uy) = C n e Dk ^ 2 gi?\m n z a ,m n z), g®(u,v) = J v (v)Y&(u) - Y v {v)J$\u), 3 = a,b, (10) 

where the normalization coefficient C n is determined below. Here J v (x), Y v (x) are the Bessel and 
Neumann functions of the order 



D 2 /A- D(D + l)( + m 2 /k 2 D , (11) 

and the z-coordinates of the branes are denoted by Zj = e k °i /ko, j = a, b. In formula (|10j) and in 
what follows for a given function F(x) we use the notation 

F (i) (x) = AjF(x) + B jX F'(x), Aj = Aj + Bjk D D/2, Bj = Bjk D , j = a, b. (12) 

Note that for conformally and minimally coupled massless scalar fields one has v = 1/2 and 
v = D/2, respectively. 

The function (|lUp satisfies the boundary condition on the brane at y = a. Imposing the 
boundary condition on the brane y = b we find that the eigenvalues m n are solutions to the 
equation 

gi ab \m n z a ,m n z b ) = {m n z a )Y^ {m n z b ) - Y^ (m n z a )jl h \m n z b ) = 0. (13) 

Equation (|13p determines the Kaluza-Klein (KK) spectrum along the transverse dimension. We 
denote by u = 7^, n = 1,2,..., the positive zeros of the function gi ab \u,uz b / z a ), arranged in the 
ascending order, 7„ jTl < j v ,n+l- The eigenvalues for m n are related to these zeros as m n = 7^/^. 
From the orthonormality condition of the radial functions, for the coefficient C n in Eq. (jlOh one 
finds 

„2 7rm n 

n ~ k D d_S ab )( vz wz ,\i ' [ ' 

du" u \ az a-, uz b)\u=m n 

Note that, as we consider the quantization in the region between the branes, a ^ y ^ b, the modes 
defined by (|10p are normalizable for all real values of v from Eq. (jlip . 

Substituting the eigenfunctions (JSj) into the corresponding mode sum and integrating over the 
angular part of the vector k, for the expectation value of the energy density on the brane at y = j 
we obtain 

ef = -In^C^zfB^ f°° dk k^ £ f) ^ , (15) 

& n=l Up,n mid" (uz a ,UZ b ) 
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where j,l = a,b, and I ^ j, and 

^ = (4 7 r)^/^r(D 1 /2)- (16) 
To regularize the divergent expression on the right of this formula we define the function 

*,-w =2^-^2,1^)1 y dfc* 1 L^,n _g, w — > ( 17 ) 

r P J{i n=l d^9iy (UZ a ,UZb)\u=m n 

where an arbitrary mass scale \i is introduced to keep the dimension of the expression. After the 
evaluation of the integral over k, this expression can be presented in the form 



where the generalized partial zeta function 



> / \ T(-a s ) ^ 2 2 m n g^\m n zi,m n Zj) Di + s-1 . . 

V ' >r" n —i Thi^v \U-Z a , UZb)\ u =m n 

The computation of the VEV of the surface energy-momentum tensor requires the analytic contin- 
uation of the function &j(s) to the value s = —1. In order to obtain this analytic continuation we 
will follow the procedure multiply used for the evaluation of the Casimir energy (see, for instance, 
|52j). On the base of this procedure it can be seen that in the strip —D\ — 1 < Res < —D\ of the 
complex plane s we have the following integral representation 

I) r(« s + 1) A, 



O/?00 = Wis) - - J - / duu{u' - \$r°n ju (uz a ,uz b ), (20) 



where 



^ T(-s/2)T{a s + l) J Xfj p Fy>(u Zj ) 



In (1201) we have defined 



i^ 6) («) li a \u) 
Q. av (u,v) = — t\ — - -, , fift^f'u, u) = -77^ — - , , (22) 

Kl a \u)Gi ab) (u,v) li b) (v)G^ b) (u,v) 

with the modified Bessel functions I v (u) and K u (u). In formula (121j) we use the notation F = K 
for j = a and F = I for j = 6. The contribution of the second term on the right of Eq. (|20|) 
is finite at s = — 1 and vanishes in the limits z a — > or z b — > oo. The first term corresponds to 
the contribution of a single brane at z = Zj when the second brane is absent. The surface energy 
density corresponding to this term is located on the surface y = a + for the brane at y = a and 
on the surface y = b — for the brane at y = b. To distinguish on which side of the brane is located 
the corresponding term we use the superscript J=L for the left side and J=R for the right side. As 
we consider the region between the branes, in formula (|2ip J=R for j = a and J=L for j = b. The 
further analytic continuation is needed for the function Qg (s) only and this is done in the next 
section. The integral representation for the partial zeta function given by Eqs. (f20|) . (|2T1) is valid 
for the case of the presence of the zero mode as well. 
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3 Surface energy density for a single brane 



In this section we consider the geometry of a single brane placed at y = j. The orbifolded version of 
this model corresponds to the higher dimensional generalization of the Randall-Sundrum 1-brane 
model with the brane location at the orbifold fixed point y = 0. The corresponding partial zeta 
function is given by Eq. (|2ip . Now in this formula J = L, R for the left and right sides of the brane, 
respectively, and F = K for J = R and F = I for J = L. In addition, the replacement — ► 
should be done with = 1 and = — 1. The integral representation (|21j) for a single brane 
partial zeta function is valid in the strip —D\ — 1 < Res < —D% and under the assumption that 
the function Fu\u) has no real zeros. For the analytic continuation to s = — 1 we employ the 
asymptotic expansions of the modified Bessel functions for large values of the argument [53] . For 
Bj 7^ from these expansions one has 

H j \u) ul ' 

where the coefficients v\ F ^{y) are combinations of the corresponding coefficients in the expansions 
for the functions F v {u) and F' u (u). Note that one has the relation = (— l)'^ , assuming 

that the coefficients in the boundary conditions are the same for both sides of the brane. For 
the nonzero modes along the internal space E we subtract and add to the integrand in (|21j) the iV 
leading terms of the corresponding asymptotic expansion and exactly integrate the asymptotic part. 
For the zero mode we first separate the integral over the interval (0, 1) and apply the described 
procedure to the integral over (l,oo). As a result, the corresponding function 



(4vr)( I) i- 1 )/ 2 ^ 



I^WQJ'OO, (24) 



is written in the form 



(47T) 



(l-D 1 )/2 n (J) z D 2 



r(-§)r(a s + i)(M*,-) s+1 



P 



+ 



duu(u 2 - \ 2 pz]) a ° 



F v (u) 
' 3 F^\u) 



N „ F >f 



E 

i=i 



i 



,i 



N 

E 

i=i 



\MX)\S 

D 1 + s - I + 1 



1 N (F,j) Dl+s _z + i 

+lr(a a + i)X: ' ^ r(- 



i=i 



r(i/2) 



a. 



FX 



(25) 



where up = ApZj + 6q\ /3 - In Eq. (p5|) we have introduced the local spectral zeta function associated 
with the massless laplacian defined on the internal subspace S: 



■2z 



(26) 



P 



where the prime on the summation sign means that the zero mode should be omitted. Both integrals 
in Eq. (|25p are finite at s = —1 for N ^ D\ — 1. For large values \p the second integral behaves 
as \p 1+s ~ N and the series over (3 in Eq. (|25p is convergent at s = —1 for N > D — 1. For these 
values iV the poles at s = — 1 are contained only in the last two terms on the right. 
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The zero mode part has a simple pole at s = — 1 presented by the summand I = D — 1 of the 
second sum in figure braces. The pole part corresponding to the nonzero modes is extracted from 
the pole structure of the local zeta function (|26|) . The latter is given by the formula 

r(*)< E (z, x)\ z=p = Cd ^ 2 -p (x) + n ( x) + • • ., (27) 

z — p 

where p is a half integer, the coefficients Cd 2 /2- p {X) are related to the Seeley-DeWitt or heat kernel 
coefficients for the corresponding non-minimal laplacian, and the dots denote the terms vanishing 
at z = p. In the way similar to that used in Ref. [23], it can be seen that the coefficients C P {X) 
are related to the corresponding coefficients C P (X, m) for the massive zeta function 



\MX)\ 

P 

by the formula 



Cds,X;rn) = ^-^-, (28) 



C P {X) = C p (X,0) - \MX)\ 2 S P ,D 2 /2- (29) 
The VEV of the energy density on a single brane is derived from 



ef = -n^k 1 j>C j <$>f{s)\ s= - l . (30) 

(J) 



By using relation (|27p . the energy density is written as a sum of pole and finite parts: 
e|j J p + . Laurent-expanding the expression on the right of Eq. (|25p near s = — 1, one finds 



for the pole part, and 



D JF,j)„ D-l 



(4vr)^/2( s + l)^ T(//2) 



~3,i 



2kBzf»c j p Dl J2\M x )f 

p 



Sox.Bj I' duu - 1 FAu) 



o 

FJu) 



I duu(u 2 -Xtff ^ [B S 



N MY 



l 



Fh J '{u) ^ 



,l 



-k°zf>C jf 3 Dl \MX)\ 2 \ E 



N ~ XF,j) 



i 2v) 



i=i 



N (F,j) 
V 



{4 ^DU-2 E f(I/2) Z f 1 \ C iD-l)/2(X) 



+ ^(/-Di)/2P0 



21n/i-V( ^ 



(32) 



for the finite part, with ip(x) being the diagamma function. In this formula the prime on the 
summation sign means that the term with I = D\ should be omitted and it is understood that 
C p (X) = for p < 0. In the pole part the second term in the square brackets comes from the zero 
mode along S and this term is cancelled by the delta term on the right of Eq. (|29p . 

The renormalization of the surface energy density can be done modifying the procedure used 
previously for the renormalization of the Casimir energy in the Randall-Sundrum model |13|. El [9] 
and in its higher-dimensional generalizations with compact internal spaces |24| [23] . The form of the 
counterterms needed for the renormalization is determined by the pole part of the surface energy 
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density given by Eq. (|31|) . For an internal manifold with no boundaries, this part has the structure 
Y^iLo *^ 2 ' a \ S \ z j I ' L) 21 . By taking into account that the intrinsic scalar curvature Rj for the brane 
& t V = j contains the factor (zj/L) 2 , we see that the pole part can be absorbed by adding to the 
brane action counterterms of the form 



[(D-l)/2) 



\h\ b i Rl r ( 33 ) 



1=0 



where the square brackets in the upper limit of summation mean the integer part of the enclosed 
expression. By taking into account that there is the freedom to perform finite renormalizations, we 
see that the renormalized surface energy density on a single brane is given by the formula 



[(D-l)/2] 

'-j ~ J ~ 



e (J,ren) = e (J) + £ cf\ Z] / L) 21 . (34) 



1=0 

The coefficients a in the finite renormalization terms are not computable within the framework 
of the model under consideration and their values should be fixed by additional renormalization 
conditions. 

The total surface energy density for a single brane at y = j is obtained by summing the 
contributions from the left and right sides: £j LR ^ = £^ + e j R ^- I n formulas (f3"Tj) . ([32]) we should 
take F = I for J = L and F = K for J = R. Now we see that, assuming the same boundary 
conditions on both sides of the brane, the coefficients C P (X, 0) enter into the sum of pole terms in 
the form 

f^/- 2i C D/2 _KX,0). (35) 



If the internal manifold contains no boundaries and D is an odd number, one has j 2 -i{X, 0) = 
and, hence, the pole parts coming from the left and right sides cancel out. For a one parameter 
internal space of size L the surface energy density on the brane at y = j is a function on the ratio 
L/zj only. Note that in the case of the AdS bulk the corresponding quantity does not depend on 
the brane position. To discuss the physics from the point of view of an observer residing on the 
brane, it is convenient to introduce rescaled coordinates 

X 'M = e -k D j x M^ M = 0, 1, 1. (36) 

With this coordinates the warp factor in the metric is equal to 1 on the brane and they are physical 
coordinates for an observer on the brane. For this observer the physical size of the subspace S is 
Lj = Le~ k °^ and the corresponding KK masses are rescaled by the warp factor: A^P = A j ge fc - Dj . 
Now we see that the surface energy density is a function on the ratio Lj /{l/ko) of the physical 
size for the internal space (for an observer residing on the brane) to the AdS curvature radius. 

As an application of the general results presented above, we can consider a simple example 
with S = S . In this case the bulk corresponds to the AdS/) + i spacetime with one compactified 
dimension X. The corresponding normalized eigenfunctions are ^^(X) = e 2m P x / L jyJTj with (3 = 
0,±1,±2, where L is the length of the compactified dimension. The surface energy density 
induced on the brane is obtained from general formulas by the replacements 

oo 

^I^POI 2 ^!^', Ap-^M D 2 = l, (37) 

P /3=0 



S 



where the prime means that the summand = should be taken with the weight 1/2. For the local 
zeta function from Eq. (|26p one has Ce(s,X) = 2L 2s - 1 Cr(2s)/(2tt) 2s , where ( R (z) is the Riemann 
zeta function. Now the only poles of the function T(z)^(z, X) are the points z = 0, 1/2. By using 
the standard formulas for the gamma function and the Riemann zeta function (see, for instance, 
[53]). it can be seen that one has Ci/ 2 (X) = —1/L, Cq(X) = 1/2^/tt, for the residues appearing in 
(|27|> and 

MX) = 2^£, jW = si p( x) = ^r(p)o,(%.), p^o, i, (38) 

for the finite parts, with 7 being the Euler constant. 



4 Two-brane geometry and induced cosmological constant 

As it has been shown in Section [21 the partial zeta function related to the surface energy density 
on the brane at y = j is presented in the form (120p , where the second term on the right is finite at 
the physical point s = — 1. For two-brane geometry the VEV of the surface energy density on the 
brane at y = j is presented as the sum 

ef=e^+Aef. (39) 

The first term on the right is the energy density induced on a single brane when the second brane 
is absent. The second term is induced by the presence of the second brane and is given by the 
formula 

. . roo 

Aef = 2C j n^\k D z j ) D B 2 ] f3 Dl ^ |#(V)| 2 / duu{u 2 - Xl) D ^ 2 - l n jv {uz a ,uz b ). (40) 

As we consider the region a ^ y ^ b, the energy desnity (I40j) is located on the surface y = a + 
for the left brane and on the surface y = b — for the right brane. Consequently, in formula (|39p 
we take J = R for j = a and J = L for j = b. The energy densities on the surfaces y = a — 
and y = b + are the same as for the corresponding single brane geometry. The expression on 
the right of Eq. (|40p is finite for all nonzero distances between the branes and is not touched by 
the renormalization procedure. For a given value of the AdS energy scale ku and one parameter 
manifold X with the length scale L, it is a function on the ratios Zb/z a and L/z a . The first ratio is 
related to the proper distance between the branes, Zb/z a = exp[/co(& — a)], and the second one is 
the ratio of the size of the internal space, measured by an observer residing on the brane at y = a, 
to the AdS curvature radius k~j^ . 

For the comparison with the case of the bulk spacetime AdS^+i when the internal space is 
absent, it is useful in addition to the VEV (1401) to consider the corresponding quantity integrated 
over the subspace E: 

Ae% = [ d^X^Aefe- ^ = e~ D ^ £ Aeg, (41) 

(s) 

where Ae^ is defined by the relation 

A<f = X>(X)| 2 A$. (42) 

& 

Comparing this integrated VEV with the corresponding formula from Ref. [32], we see that the 
contribution of the zero KK mode (A« = 0) in Eq. (|4ip differs from the VEV of the energy density 



9 



in the bulk AdSo 1+ i by the order of the modified Bessel functions: for the latter case v — > v\ with 
v\ defined by Eq. (jlip with the replacement D — > D\. 

Now we turn to the investigation of the part (|4U|) in the surface energy density in asymptotic 
regions of the parameters. For large values of AdS radius compared with the interbrane distance, 
— a) "C 1, the main contribution to the integral on the right of Eq. (|40p comes from large 
values of uz a ~ [kn(b — a)]" 1 . Assuming that B a /(b — a) and m(b — a) are fixed, we see that 
the order of the modified Bessel functions is large. Replacing these functions by their uniform 
asymptotic expansions for large values of the order [53], it can be seen that to the leading order 
the corresponding surface energy on the branes in the bulk geometry ij(- Dl_1 ' 1 ) x S is obtained. 

For large KK masses along S, z a \p 3> 1, A^ 3> 1, we can replace the modified Bessel functions 
by the corresponding asymptotic expansions for large values of the argument. For the contribution 
of a given KK mode to the leading order this gives 



where 



I 42 D2,,2, y 2\-l/_.2 _ \ 

A £ .^4n k D Zj B 3 C 3 p Dl J^duu Ca{uzMuZb)eMzb - Za) _ x > 



£1 
2 



Cj ( u ) = -± rr^, j = a, b. (44) 

JV ' Aj+n^Bju v ' 

If in addition one has the condition \p{z b — z a ) 3> 1, the dominant contribution into the w-integral 
comes from the lower limit and we have the formula 

M 2n^BlCj k°z» +l X ^/^ e -2X,{ Zb - Za ) 

>*> ^-(A^ JJ B,) 2 (4vr) B i/ 2 Ca (A^ a ) Cb (A^ & )(^-z a )^/ 2 ' 1 ; 

In particular, for sufficiently small length scale of the internal space this formula is valid for all 
nonzero KK masses and the main contribution to the surface densities comes from the zero KK 
mode. In the opposite limit of large internal space, to the leading order we obtain the corresponding 
result for parallel branes in AdSz) + i bulk [32] . 

For small interbrane distances, kr>{b — a) <C 1, which is equivalent to z b /z a -1<1, the main 
contribution into the integral in Eq. (|40p comes from large values u and to the leading order we 



obtain formula (|43H . If in addition one has Xp{z b — z a ) <C 1 or equivalently Ai (6 - o)<l, we can 
put in this formula \p = 0. Assuming (b — a) <C \Bj/Aj\, for Bj 7^ to the leading order one finds 

where Gj = 1 for \B a /A a \, \B b /A b \ > ko(b - a), and <Tj = 2 2 " D - 1 for \Bj/Aj\ > k D (b - a) and 
B[/A[ = 0, with I = b for j = a and I = a for j = b. We see that for small interbrane distances 
the sign of the induced surface energy density is determined by the coefficient Cj and this sign is 
different for two cases of aj. 

Now we consider the limit Xpz b 3> 1 assuming that \pz a < 1. Using the asymptotic formulas 
for the Bessel modified functions containing in the argument z b , for the contribution of a given 
nonzero KK mode we find the following results 



{s) _^ k»z°BlC a (\p/z b ) D ^e-^*» 
* ±7r C b (\pZ b )K y v {\pz a ) 



k D 7 D + 1 B' 2 r l . \? l/2+1 e- 2X Zb T {a) (\o7 \ 
W ~ 2^-H D ^-^ / HA b + \pz b B b yKi a \\ f3 z a ) [ ' 
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This limit corresponds to the interbrane distances much larger compared with the AdS curvature 
radius and with the inverse KK masses measured by an observer on the left brane: (b — a) S> 
1/ko, 1/X^a\ For a single parameter manifold £ with length scale L and (b — a) 3> L a these 
conditions are satisfied for all nonzero KK modes. 

In the limit z a X/3 <C 1 for fixed z b A/3, by using the asymptotic formulas for the modified Bessel 
functions for small values of the argument and assuming \A a \ ^ \B a \i>, one finds 

A (s) _ k D z a +2U BlC a (j Dl [°° j 2u+l( 2 \2\£>i/2-l jj^ \uz b ) 



*® ~ (50) 



where we have introduced the notation 

(6) _ 4BZC b [3 Dl r» du J^ (u 2 _ X 2 z 2 )Dl /2-l m 

For \A a \ = \B a \u we should take into account the next terms in the corresponding expansions of 
the modified Bessel functions. The integral in Eq. (|49p is negative for small values of the ratio 
A b /B b and is positive for large values of this ratio. As it follows from Eq. (|50l) . for large interbtane 
separations the sign of the quantity Ae[ s ^ is determined by the combination {B 2 u 2 — A 2 )C b of the 
coefficients in the boundary conditions. In the limit under consideration the KK masses measured 
by an observer on the brane at y = a are much less than the AdS energy scale, C i]), and 
the interbrane distance is much larger than the AdS curvature radius. In particular, substituting 
\p = 0, from these formulas we obtain the asymptotic behavior for the contribution of the zero 
mode to the surface energy density induced by the second brane in the limit z a /z b -C 1. Now 
combining the corresponding result with formulas (|47p . (|48p . we see that under the conditions 
(6— a) 3> l/kD,L a , and L a ko ^ 1, the contribution of the nonzero KK modes along £ is suppressed 
with respect to the contribution of the zero mode by the factor (z b / ' z a ) 2u+Dl ^ 2+S i exp(— 2\pz b ). 

From the analysis given above it follows that in the limit when the right brane tends to the AdS 
horizon, z b — ► oo, the energy density Ae^g vanishes as e~ 2X ^ Zb / z^ 1 for the nonzero KK mode 

along S and as z b ' Dl ^ 2l/ for the zero mode. The energy density on the right brane, Ae^g, vanishes 
as Z ^ 2+Dl l 2+1 e -2\pz h £ Qr nonzero KK mode and behaves like z® 2 ~ 2v for the zero mode. In 
the limit when the left brane tends to the AdS boundary, z a — > 0, the contribution of a given KK 
mode vanishes as z® +2u for Aejjfg and as z 2u for Ae^g. For small values of the AdS curvature 
radius corresponding to strong gravitational fields, assuming \pz a 3> 1 and — z a ) 3> 1, we 

can estimate the contribution of the nonzero KK modes to the induced energy densities by formula 
(I45p . In particular, for the case of a single parameter internal space with the length scale L, under 
the assumed conditions the length scale of the internal space measured by an observer on the 
brane at y = a is much smaller compared to the AdS curvature radius, L a -C k^ 1 . If L a > 
one has \gz a < 1 and to estimate the contribution of the induced surface densities we can use 
formulas (j47]) and (I48p . and the suppression is stronger compared with the previous case. For 
the zero KK mode, under the condition kz>(b — a) 3> 1 we have z a /z b <C 1 and to the leading 
order the corresponding energy densities are described by relations (I49p and (|50p . From these 
formulae it follows that the induced energy densities integrated over the internal space behave as 
kp 1+1 exjp[(Di5j + 2v)k£)(a — b)] for the brane at y = j and are exponentially suppressed. 

Introducing the rescaled coordinates defined by Eq. (|36p . after the Kaluza-Klein reduction of the 
higher dimensional Hilbert action, by the way similar to that in the Randall-Sundrum braneworld, 
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it can be seen that effective Di-dimensional Newton's constant Gd 1 j measured by an observer on 
the brane at y = j and the fundamental {D + l)-dimensional Newton's constant Gd+i are related 
by the formula 

<7n (D -2)k D G D+1 ( D -2)k n (b-i) y -D 2 k D j f jD 2 y r~ /j-gv 

3 ~ V^j [ e (D-2)k D (b-a) _ x] e ' V *3 ~ e J d A V7, 

where Vjy is the volume of the internal space measured by the same observer. In the orbifolded 
version of the model an additional factor 2 appears in the denominator of the expression on 
the right. Formula (|52p explicitly shows two possibilities for the hierarchy generation by the 
redshift and large volume effects. Note that the ratio of the Newton constants on the branes, 
Goib/GD 1 a = e^ 2 ~ Dl ^ kD ^ b ~ a \ is the same as in the model without an internal space. For large 
interbrane distances one has Gz> ia ~ koGo+i/Vsa, Gd^ ~ k^GD+i^ 2 ^ ^ ^^ /Vsb, and the 
gravitational interactions on the brane y = b are exponentially suppressed. This feature is used in 
the Randall- Sundrum model to address the hierarchy problem. As we will see below this mechanism 
also allows to obtain a naturally small cosmological constant generated by the vacuum quantum 
fluctuations (for the discussion of the cosmological constant problem within the framework of 
braneworld models see references given in |32j). 

As we have already mentioned, surface energy density (|40p corresponds to the gravitational 
source of the cosmological constant type induced on the brane at y = j by the presence of the second 
brane. For an observer living on the brane at y = j the corresponding effective Di-dimensional 
cosmological constant is determined by the relation 

A Dlj = SttGd^Ae^ = 8nM 2 Di ^A8^, (53) 

where Md 1 j is the Di-dimensional effective Planck mass scale for the same observer and • 
is defined by Eq. (|4ip . Denoting by M^+i the fundamental (D + l)-dimensional Planck mass, 
G D+1 = M^g, from Eq. © 

one has the following relation 

M^\ Dl - 2 _ (z b /z a ) D - 2 - 1 D2+1 

M D+l J ~ (D - 2)(z b / Zj )D-2 kD M d + i > ^ 

for the ratio of the effective and fundamental Planck scales. By using the asypmtotic relations given 
above, for large interbrane distances one obtains the following estimate for the ratio of the induced 
cosmological constant (j53[) to the corresponding Planck scale quantity in the brane universe: 



For the model without an internal space this ratio is of the same order of magnitude for both 
branes. 

In the higher dimensional version of the Randall- Sundrum braneworld the brane at z = z b 
corresponds to the visible brane. For large interbrane distances, by taking into account Eq. (|50p . for 
the ratio of the induced cosmological constant (|53p to the Planck scale quantity in the corresponding 
brane universe one obtains 

where the function ffty is defined by Eq. ([5Tj) . 



-2 
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Using relation (|54|) with j = b, we can express the corresponding interbrane distance in terms 
of the ratio of the Planck scales 



M Di+X v ^b 
"/) • 1 



k D 



Mp+l 
M Dlb 



1 

D-2 



Substituting this into Eq. (|56p . for the ratio of the cosmological constant on the brane at j 
the corresponding Planck scale quantity one finds 



h b 



1 



caiy) \M D+ iJ 



Di-0 



D L +u(D 1 -2) 



f(b) 
Jv/3 ' 



(57) 



b to 



(58) 



with v = 2v/(D — 2). The higher dimensional Planck mass Md+i and AdS inverse radius kr> 
are two fundamental energy scales in the theory which in the Randall-Sundrum model are usually 
assumed to be of the same order, kn ~ Md+i- In this case one obtains the induced cosmological 
constant which is exponentially suppressed compared with the corresponding Planck scale quantity 
on the visible brane. In the model with D\ = 4, hp ~ Md+i ~ 1 TeV, Md-,6 = Mp\ ~ 10 16 
TeV, assuming that the compactification scale on the visible brane is close to the fundamental 
Planck scale, Vs6-^d+i ~ ^' ^ or ^ ne ra ti° °f the induced cosmological constant to the Planck scale 
quantity on the visible brane we find the estimate h b ~ 10~ 32 ( 2 + U ) . From (|57p one has ko{b — a) ~ 
74/ (D2 + 2) and the corresponding interbrane distances generating the required hierarchy between 
the electroweak and Planck scales are smaller than those for the model without an internal space. In 
the model proposed in Ref. [23], a separation between the fundamental Planck scale and curvature 
scale is assumed: ko ~ Mo+iZa/ Zb ~ 1 TeV. Under the assumption Vs&Mm-i ~ 1> ^ n this model 
we have h b ~ i - 64 [ 1 +^( D + 1 )l and k D {b - a) w 74/(D 2 + 3). 



5 Conclusion 

We have investigated the expectation value of the surface energy-momentum tensor induced by the 
vacuum fluctuations of a bulk scalar field with an arbitrary curvature coupling parameter satisfying 
Robin boundary conditions on two parallel branes in background spacetime AdS£> 1+ i x S with a 
warped internal space S. Vacuum stresses on the brane are the same for both subspaces and the 
energy-momentum tensor on the brane corresponds to the source of the cosmological constant type 
in the brane universe. It is remarkable that the latter property is valid also for the more general 
model with the metric instead of rj^ in line element ([2]). As an regularization procedure for 
the surface energy density we employ the zeta function technique. The corresponding zeta function 
is presented as the sum of single brane and second brane induced parts. The latter is finite at the 
physical point and the further analytical continuation is necessary for the first term only. As the 
first step we subtract and add to the integrand the leading terms of the corresponding asymptotic 
expansion for large values of the argument and explicitly integrate the asymptotic part. Further, for 
the regularization of the sum over the modes along the internal space we use the local zeta function 
related to these modes. By making use of the formula for the pole structure of this function, we 
have presented the energy density on a single brane as the sum of the pole and finite parts. The pole 
parts in the surface energy density are absorbed by adding to the brane action the counterterms 
having the structure given by Eq. (|33p . The renormalized energy density on the corresponding 
surface of a single brane is determined by formula (|34p , where the second term on the right presents 
the finite renormalization part. The coefficients in this part cannot be determined within the model 
under consideration and their values should be fixed by additional renormalization conditions which 
relate them to observables. 
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Unlike to the single brane part, the surface energy density induced by the presence of the 
second brane contains no renormalization ambiguities and is investigated in Section [H This part 
is given by formula (|40|) . We have investigated the induced energy density in various asymptotic 
regions for the parameters of the model. In the model under discussion the hierarchy between the 
fundamental Planck scale and the effective Planck scale in the brane universe is generated by the 
combination of redshift and large volume effects. The corresponding effective Newton's constant 
on the brane at y = j is related to the higher-dimensional fundamental Newton's constant by 
formula f)52 j) and for large interbrane separations is exponentially small on the brane y = b. We 
show that this mechanism also allows obtaining a naturally small cosmological constant generated 
by the vacuum quantum fluctuations of a bulk scalar. For large interbrane distances the ratio 
of the induced cosmological constant to the the corresponding Planck scale quantity in the brane 
universe is estimated by formula (|55f) and is exponentially small. For the visible brane in the 
higher dimensional generalization of the Randall-Sundrum two-brane model, this ratio is given in 
terms of the effective and fundamental Planck masses by Eq. (|58|) , We have considered two classes 
of models with the compactification scale on the visible brane close to the fundamental Planck 
scale. For the first one the higher dimensional Planck mass and the AdS inverse radius are of the 
same order and in the second one a separation between these scales is assumed. In both cases the 
corresponding interbrane distances generating the hierarchy between the electroweak and Planck 
scales are smaller than those for the model without an internal space and the required suppression 
of the cosmological constant is obtained without fine tuning. 
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